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By using Majorana’s stellar representation, we give a clear geometrical interpretation of the topo¬ 
logical phases of inversion-symmetric polymerized models by mapping the Bloch states of multi-band 
systems to Majorana stars on the Bloch sphere. While trajectories of Majorana stars of a filled Bloch 
band exhibit quite different geometrical structures for topologically trivial and nontrivial phases, 
we further demonstrate that these structures are uniquely determined by distributions of Majo¬ 
rana stars of two high-symmetrical momentum states, which have different parities for topologically 
different states. 
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I. INTRODUCTION 

Recently, Majorana’s stellar representation (MSR) has 
gained revived interest as it offers an intuitive geomet¬ 
rical approach to understanding quantum systems with 
multiple components by mapping states in the higher di¬ 
mensional space as points on a Bloch sphere^'— . In gen¬ 
eral, a quantum state of spin-1/2 systems or equivalently 
two-level systems can be represented as a point on a unit 
sphere^, and the evolution of the quantum state corre¬ 
sponds to a particular trajectory on the Bloch sphere. 
By using MSR, which represents a quantum pure state 
of spin-J systems in terms of a symmetrized state of 2J 
spin-1/2 systems, one can generalize this geometric ap¬ 
proach to large spin systems or multi-level systems. In 
terms of MSR, the evolution of a spin-J state can be 
intuitively understood by trajectories of 2J points on 
the two-dimensional (2D) Bloch sphere, with these 2J 
points generally coined as Majorana stars (MSs). This 
approach naturally provides an intuitive way to study 
high spin systems from geometrical perspectives, which 
has made the MSR a useful tool in many different fields, 
e.g., classification of entanglement in symmetric quan¬ 
tum states^-—, analyzing the spectrum of the Lipkin- 
Meshkov-Click modeU^, studying Bose condensate with 
high spinsi^“— , and calculating geometrical phases of 
large-spin systems^i^. 

As much of the search for the application of MSR has 
focused on high-spin systems, it is interesting to apply 
this approach to study the multi-band topological sys¬ 
tems. It is well known that a topological insulator distin¬ 
guishes a trivial band insulator by its nontrivial topolog¬ 
ical energy band, which has different geometrical prop¬ 
erty from a trivial bandit. For a two-band system, e.g., 
the Su-Schrieffer-Heeger (SSH) modeU^, one can map the 
Bloch state into a 2D Bloch sphere, and the geometrical 
meaning of topologically different phases can be unveiled 


by their distinct trajectories^Si^. As a paradigmatic 
topological model2^, the SSH model supports either topo¬ 
logically trivial or nontrivial phase, characterized by the 
quantized Berry phase 0 or which is experimen¬ 

tally measurable as demonstrated in the recent cold atom 
experiments^. In this work, we shall apply MSR to study 
one-dimensional (ID) topologically nontrivial polymer¬ 
ized systems, which can be viewed as a multi-band gen¬ 
eralization of the celebrated dimerized SSH modeUSiSi 
and were also proposed to be realizable in optical su¬ 
perlattice systemsSSiSS. The application of MSR enables 
us to investigate the geometrical property of muti-band 
topological systems by mapping the multi-level quantum 
states to MSs on the Bloch sphere. Consequently, a filled 
Bloch band forms specific trajectory of MSs on the Bloch 
sphere, which exhibits very different geometrical struc¬ 
ture for topological or trivial phase. Furthermore, we 
unveil that the distinct geometrical structures of trajecto¬ 
ries for topologically different states are closely related to 
their parities and determined by the distribution of MSs 
at two high-symmetry momentum points. Our study pro¬ 
vides an intuitive way to distinguish topologically differ¬ 
ent phases of multi-band systems and unveils the intrinsic 
relation between the band topology and parity from the 
geometrical point of view, which shall deepen our under¬ 
standing of their geometrical properties. 


II. ID POLYMERIZED MODEL WITH 
INVERSION SYMMETRY 

We consider a ID superlattice with period T (T > 2) 
described by the Hamiltonian: 

T 

H = ^ ^ ^ ^ ^ h.C.), (1) 

n a=l 


2 


where CT+i,n = and we require the system having 

inversion symmetry by enforcing 


ta — ^T—c 


( 2 ) 


For the case with T = 2, it reduces to the SSH modeU^. 
Models with T > 3 can be viewed as generalizations 
of the dimerized SSH model^i, and for convenience, we 
also refer them as polymerized models. As the unit 
cell consists of T sites, the spectrum of the system 
splits into T bands. By taking a Fourier transforma¬ 
tion Ca^n = 1/V^ with a = 1,2,-•• ,T, the 

Hamiltonian in the momentum space can be written as 
H = T,k with ipk = (cife, C2fc, • • • , CTkV and 
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not have topologically nontrivial properties. Neverthe¬ 
less, when the system has inversion symmetry, it may 
support topologically nontrivial phase characterized by 
the quantized nontrivial Berry phase or the emergence of 
doubly degenerate edge states for the open chain^S. For 
convenience, we shall focus our study on a concrete case 
with T = 3 in the present work, and also show our results 
can be directly generalized to cases with larger periods, 
e.g., the case with T = 4. 

For the T = 3 superlattice model or trimerized model 
with inversion symmetry, the hopping amplitudes can be 
parameterized as 


ti =t2 = t{l - 6), ts = t{l + (5), 


( 4 ) 


by Pij = The system also has time reversal 

symmetry, which leads to Th{k)T~^ = h{—k), here the 
time reversal operator T is just the complex-conjugation 
operator K. As we shall see in the following context, both 
the inversion and time reversal symmetries give some re¬ 
strictions on the MSR of Bloch states, which plays an 
important role in determining trajectories of MSs of topo¬ 
logically different states. 


III. MAJORANA REPRESENTATION AND 
BERRY PHASE 

For a T-band system, the Bloch state can be expressed 
as |4>(/c)) = ^^=iCa{k)\a)k> To represent this multi¬ 
level state by MSs, it is convenient to map the state to 
a spin-J state |4>(/c)) = ^ ~ 

(T — l)/2. There is a one-to-one correspondence between 
parameters Ca{k) and Cm{k) by taking m = a — 1 — 
J. According to Schwinger boson representation theory^, 
the angular momentum operators can be described by 
creation and annihilation operators of two mode bosons, 
d+, d, and 6 , and the state |J, m) can be expanded 
by |J,m) = [(J + m)!(J — 

where jp) is defined by d+jp) = | t) and 6 +|p) = | |). 
With the help of Schwinger representation, the state 4>(/c) 
can be factorized as 

2 J 

=W (5) 

^ i=i 

where Nj is the normalization coefficient. If we denote 
Xj = tan , then the factorization parameters Oj and 

(pj can be determined by the roots of the following poly¬ 
nomial equatioi>i 


where t = 1 is taken as the unit of energy and |d| < 1 is 
set. The spectrum is split into 3 bands and gap between 
bands is always open for any nonzero 5. Considering the 
state with the lowest band being fully filled, we find that 
the state is topologically different for d > 0 and d < 0 . 
For the open chain with the length L = 3A, there appear 
degenerate edge states at both ends for d > 0 but none for 
d < 0. This result suggests that there exists a topological 
phase transition by varying d with the transition point 
at d = 0. For the periodic chain, the topological phase 
transition can be characterized by the change of the Berry 
phase of the system, i.e., the Berry phase 7 = tt in the 
topological phase and 7 = 0 in the trivial phase. Here the 
Berry phase across the Brillouin zone (BZ), also known 
as Zak phased!, is defined as 7 = i §\^{k)\-^\^{k))dk 
with 4>(/c) denoting the occupied Bloch states. 

The existence of topological states in our polymerized 
models is protected by inversion symmetr y^^i^^ . In the 
momentum space, the inversion symmetry means that 
Ph{k)P~^ = h{—k), where the inversion operator P is an 
anti-diagonal matrix with the the matrix element given 


(-ircj-j(fc) 2j-j _ n 


(6) 


From Eq.(j5j), it is obvious that |4>(/c)) can be viewed 
as the product of 2 J spin-^ states with \uj) = 

^cos ^, sin ^ , and a given T-band Bloch state can 

be described by T — 1 MSs on the Bloch sphere. 

Considering the time reversal symmetry and inversion 
symmetry of the inversion-invariant polymerized model, 
we can directly get the eigenstate |4>) fulfilling the fol¬ 
lowing relations: 


\^*{k)) = m-k)), ( 7 ) 

\<^*{k)) = p\m)- ( 8 ) 

From Eq.Q, we get Cj{k) = Cj{—k). So if x{k) = 
tan is the solution of Eq.(j 6 ]), then x*(/c) must be the 
solution of Cj-j{-k)/y/{2J - = 

0. In other words, for any star Uj{k) of {0j{k)^(j)j{k))^ 
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there must exist a star ui{—k) satisfied Oj{k) = Oi{—k), 
(t^j{k) = —(j)i{—k). So the whole trajectory of MSs for a 
filled Bloch band is symmetric about the meridian of the 
Bloch sphere. Similarly, from Eq.([ 8 ]), we have Cj-j{k) = 
Cj_j{k), j = 0 , 1 , • • • , n. So if X = tan is a solution 

of Eq.(j 6 ]), x' = = tan is also a solution. That 

is to say, for a fixed momentum k, any star Uj of {0j^(j)j) 
corresponds to a star ui of {0i = ir — Oj, (pi = 0 j), unless 
Uj lies on the equator. So the whole trajectory of MSs 
for a filled band is always symmetric about the equator. 



FIG. 1: Trajectories of MSs for the lowest Bloch band of the 
trimerized model with inversion symmetry. The blue and red 
curves represent different MSs. (a) topologically trivial state 
with 6 = —0.2, the Berry phase is 71 = 0.404, 72 = —0.404, 
7 c = 0. (b) topologically nontrivial state with S = 0 . 2 , the 
Berry phase is 71 = —2.216, 72 = —0.925, jc = 0. (c) and 
(d) show the change of <p of MSs for the lowest Bloch band of 
trimerized model as k goes from — tt to tt, corresponding to 
(a) and (b), respectively. 


In terms of the MSR, the Berry phase for the T-band 
polymerized model can be represented as a sum of two 
parts, 7 = 70 + 7c, where 

70 = E “ E ^ f , 

j 3 


and the phase 



Ui X Uj • d{uj — Ui) dNj 


represents correlations between any two MSs^, where Ui 
is a function of k and the integral is carried out along 


the evolution paths of MSs as the momentum goes over 
the BZ. By using symmetries of MSs, we can prove that 
7 c = 0 as the summation terms of correlations cancel out 
each other (see the appendix for details). Consequently, 
the Berry phase of the occupied Bloch band is simplified 
to the sum of Berry phases of MSs on the Bloch sphere, 
generated by the evolution of each star. Since the MSs 
either appear in pairs, in terms of Uj and ui{0i = ir — 
Oj^pi — (j)j) for Oj y or locate in the equator, we can 
further simplify the expression of Berry phase by using 
'^j cos Oj = 0. Einally, the Berry phase of an occupied 
Bloch band is simplified to 

which is the sum of integrals of each MS along its pro¬ 
jecting trace on the equator. 

Eor the trimerized model parameterized by Eq.(|4]), the 
eigenstate |^(/c)) can be described by two MSs. As k 
goes across the BZ, trajectories of the lowest Bloch band 
for cases with < 0 and ^ > 0 are shown in Eig(T]^a) and 
(b), respectively. It is obvious that trajectories for topo¬ 
logically different states display distinct geometric struc¬ 
tures. Here, in the topologically trivial phase, the Berry 
phase of each star cancels out. However, in the topo¬ 
logically nontrivial phase, the trajectories of two stars 
splice together to cover the whole equator, which gives 
7 = TT. To see this more clearly, we show the change 
of as a a function of momentum in Eig(T]^c) and (d), 
corresponding to EiglTJa) and (b), respectively. Eor the 
topologically trivial case, as k goes from — tt to tt, the in¬ 
tegral over (j) for each star cancels out each other, which 
generates a zero Zak phase. 

Next we show that different topological structures of 
trajectories of MSs are closely related to properties of 
MSs at /c = 0 and k = tt. At these high-symmetry 
points denoted by kg, which fulfils ks up to a re¬ 

ciprocal lattice vector, we have Th{ks)T~^ = h{ks) and 
Ph{ks)P~^ = h{ks). These symmetries suggest that the 
distribution of MSs for the high-symmetry state ^{ks) 
must be mirror-symmetrical about both the meridian and 
equator. The state ^{ks) also has a certain parity, i.e., 
P\^{ks)) = ^{ks)\^{ks)) with ^{ks) = ± 1 , and the parity 
^ of a filled Bloch band is given by the product of ^(0) 
and ^(tt). Using the expression of P, we get 

Cj-j{ks)=i{ks)Cj-j{ks), ( 10 ) 

= 1 , ( 11 ) 

where all coefficients Cj{ks) are real due to the time- 
reversal symmetry. Substituting m into the Eq.( 6 ), we 
find that if Xj is the solution of the equation, then x*, 
and (^)* are also the solutions, i.e., the Majorana 
stars for the high-symmetry state are mirror-symmetrical 
about both the meridian and equator, as schematically 
shown in Eig. [2j If the four stars are degenerated 

to 2 stars on the equator; if pj = 0 or tt, the stars are 
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degenerated to two stars on the meridian; and if Oj = ^ 
and (j)j = 0 or TT, all the stars are degenerated to a single 
star. For the odd-parity state, ^{kg) = —1, and we get 
(pj = TT from Eq. ([TT]) . While for the even-parity 
state, ^{ks) = 1, and we get (j)j = 0. 


= 0 



FIG. 2: Schematic plot for the distribution of Majorana 
stars of high-symmetry momentum state. At the high- 
symmetry point, every four stars form a group which is 
mirror-symmetrical about both the meridian and equator. 

For the trimerized model with T = 3, there are only 
two Majorana stars corresponding to a high-symmetry 
momentum state. For the odd-parity state, the stars 
must be located at 'Ui(|^,0) and 'U2(f,7r) in order to 
fulfil the mirror-symmetry conditions and the condition 
of + 02 = tt simultaneously. While for the even- 
parity state, the stars can be located at 'Ui(|^,0) and 
'^ 2 (f 5 “ 0 ) or at uiiO^Q) and 'U2(7r —^,0), or at 'Ui(^,7r) 
and 142 (tt — To see it clearly, we display the distri¬ 

bution of MSs at the high-symmetry points for both the 
topologically trivial (in FiglSl^a)) and nontrivial phases 
(in FiglUb)). While the parity of 4>(7r) is always even 
with ^(tt) = 1 for both phases, we have ^(0) = 1 for the 
trivial phase, and ^(0) = —1 for the topological phase. 

As the topologically trivial and nontrivial phases have 
different parities, next we unveil the Berry phase of the 
corresponding Bloch band can be uniquely determined 
by distributions of MSs at two high-symmetry points. 
To calculate the Berry phase, it is convenient to project 
the trajectory of MSs on the Bloch sphere to the equa¬ 
tor, according to Eq. Consequently the MSs sym¬ 
metrically located above and below the equator merge 
to the same point. Then, we divide the BZ into two 
parts ( 0 , 7 r) and (7r,27r), and we just need to consider 
the interval ( 0 , 7 r) as the integral of 0 over the interval 
(tt, 27r) (equivalently (—tt, 0 )) gives the same contribution 
to the Berry phase as over the interval ( 0 , 7 r), according 
to the symmetry analysis. Hence, the Zak phase 7 is 
just the double of the Berry phase as k integrates from 
0 to TT. In the topologically trivial phase, when /c = 0, 
both MSs are located at “Oi” on the projected equator 
shown in EigHJ^a). As k travels from 0 to tt, the stars 




FIG. 3: MSR for the high-symmetrical momentum state of 
the trimerized model with (a) S = —0.2 and (b) S = 0.2, 
respectively. Blue points and red stars correspond to the state 
at /c = 0 and /c = tt, respectively. While the parity for the 
point of /c = TT is always even, the parity for the point of /c = 0 
can be even for the trivial phase (a) or odd for the topological 
phase (b). 



FIG. 4: Schematic plot for the change of polar angle of MSs 
when k goes from 0 to tt, (a) for topologically trivial phase, (b) 
and (c) for topologically nontrivial phase. Here Oj correspond 
to MSs at /c = 0 and Hj to MSs at /c = tt. 


move to Hi and n 2 , generating changes of polar angles 
A 01 = 2 m 7 r + ai and A 02 = 2 n 7 r — o^i, respectively, 
where m^n can be any integer. From Eq.(|9]), we get 
the Zak phase 7 = —(A0i + A 02 ) = —2{m + n) 7 r, i.e., 
7 mod (27r) = 0. In the topologically nontrivial phase, 
MSs corresponding to the /c = 0 state are located at 
“Oi” and “ 02”5 and there are two possibilities for the 
evolution of state from /c = 0 to /c = tt, as shown in 
FiglH^b) and (c). For the case of FigHJ^b), Oi moves to 
Hi and O 2 to n 2 as k goes from 0 to tt, with correspond¬ 
ing changes of polar angles given by A0i = 2 m 7 r + 0^2 
and A 02 = 2 n 7 r -h 02 - Since 0^2 + 02 = tt, we get the Zak 
phase 7 = —2{m + n) 7 r — tt, i.e., 7 mod (27r) = tt. For 
the case of FiglH^c), we can get the same conclusion by 
following similar discussions. 

Our results can be directly generalized to general poly¬ 
merized models with inversion symmetry. Concretely, 
we consider a tetramerized model with ti = ts = t, 
t 2 = t{l — S) and ^4 = t(l +J), and take t = 1 and |J| < 1 . 
Similar to the case of T = 3, there exists a topological 
phase transition occurring at (5 = 0, with (5 < 0 corre¬ 
sponding to the trivial phase and > 0 the topologically 
nontrivial phase. The analysis of geometric meaning is 
analogous to the trimerized model despite a little more 
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complicated. 



FIG. 5: MSR for the topologically trivial and nontrivial state 
of the tetramerized model with (a) S = —0.2 and (b) S = 
0.2, respectively. The blue, red and green curves represent 
different MSs. Blue points and red stars represent MSs of 
states at /c = 0 and /c = tt, respectively, (c) and (d) display 
the change of 0 of MSs as k goes from — tt to tt, corresponding 
to (a) and (b), respectively. 

For the tetramerized model with T = 4, each state in 
the momentum space is represented by three MSs, and 
their trajectories over the BZ are shown in FigISfa) and 
(b) for topologically trivial and nontrivial phases, respec¬ 
tively. MSs in two high-symmetry points are also marked 
on the trajectories. In the topologically trivial phase, 
all trajectories of MSs go back and forth, whereas in 
the topologically nontrivial phase, the trajectories form 
a close circle and cover the equator completely. Conse¬ 
quently, Zak phases for the topologically trivial and non¬ 
trivial states are 0 and tt, corresponding to Figlbfc) and 
(d), respectively. Similarly, different MS distributions at 
k = TT for topologically trivial and nontrivial states indi¬ 
cate that they have different parities. 


IV. SUMMARY 

In summary, we have unveiled the geometrical proper¬ 
ties of topological phases of ID inversion-invariant multi¬ 
band systems from trajectories of MSs, which exhibit 
different topological structures for topologically different 
phases. By utilizing the symmetric properties of MSs, we 
found that the Zak phase of a filled band can be repre¬ 
sented as the summation of integral of each MS along its 
projecting trace on the equator, and takes 0 and tt for 
topologically trivial and nontrivial phases, respectively. 
We further demonstrated that the topological structure 


of the trajectory of MSs is closely related to the parity of 
the system, which is determined by properties of Bloch 
states at two high-symmetry points. 
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Appendix A: The calculation of correlation phase 


The correlation phase 7c in the main text can be rep¬ 
resented as 7 c = / d^c with 



Ui X Uj • d{uj — Ui) dNj 


(Al) 


where the normalization coefficient Nj is a symmetric 
function for permuting any two stars Ui and Uj , and dij = 
1 — Ui^ Uj is the distance. It is clear that 7c is a sum of 
correlations between any two stars. Due to the existence 
of time reversal symmetry and inversion symmetry, the 
Majorana stars for a Bloch state with fixed momentum k 
distribute either on the equator or symmetrically about 
the equator. So the distribution of Majorana stars for 
a Bloch state can be schematically displayed in Fig. [6l 
where Ai are the stars above the equator, Bi are the 
corresponding ones below the equator symmetrical to Ai , 
and Ci are on the equator. 
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FIG. 6: Schematic plot for the distribution of all the Majorana 
stars on the 0 — (j) plane for a Bloch state with fixed k. The 
stars are classified into three classes denoted by Ai, Bi and 
Ci, respectively. 
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Firstly, we calculate the correlation between Ai and 
Aj and the corresponding correlation between Bi and 
Bj, which are given by 


dlcAiAj 

sin OAi cos (l)Ai sin 6^^. sin (pAi cos 6^^. 
sin Oaj cos 4>Aj sin Oa^ sin 4>Aj cos Oa^ 

ail 0112 <^13 


dk 

T 


dNj 


NjddA.A, 


Since 6^^. = tt — Osi and (pAi = <pBi^ we have = 


dsi 


N'iddA 


and <P'a, = 


0^., where 0 = dO/dk and p = d<p/dk. Substituting 
them into (ixii and dssi we find that d^cAiAj -^d^cBiBj = 
0. Similarly, we can derive djcAiBj + djcBiAj = 0, 
djcAiCi + djcBiCi = 0. 

For the correlation between Ai and Bi, we know that 
(A2)^^i = tt — OBi = 0 and pAi = <pBi = 0, so 


and 


djcBiBj 


& 

T 


sin 0Bi cos (j)Bi sin 0Bi sin pBi cos 0Bi 
sin Obj cos (pBj sin Ob^ sin pBj cos Ob^ 
a2l a22 <^23 


dNj 

NjddB^B, ’ 


d'JcAiBi 

sin 6* cos ^ sin 0 sin ^ cos 6 
sin 9 cos (f> sin 9 sin (f> — cos 9 
0 0 26»'sin6> 

= 0 . 


Y 


dNj 

NjddA.B, 


(A3) 

where 

an = (cos6*A,- cos4 >a^9'a. - sin6'A^ sin^A^X) “ 0 

ai2 = (cos6»a,. sin</>A3 6*A, + sin^yi, cos</>AY/ij) “ 

ai3 = -sin^A, A- “ ^ 

a 2 i = {cos 9 b j cos 4 )b^9b. - sin 6»Bj sinc/iB^X) “ (•? ^ 

a 22 = (cos 6'b, sin ^b, 9b. + sin 6 'b, cos (/)Bj ) - {j ^ i), 

a23 = -sin^B. ^B, - (j ^ 


At last, the correlation between Ci and Ck is always 0, 
because they are on the same plane so that ui x Uk^d{ui — 
Uk) vanishes apparently. Hence, we have 

d'ye = {d-ycAiA^ + d'ycBiBj) + {d'ycAiBj + d'ycBtA^) + 
{d'ycAiCi + d'ycBiCi) + d'ycAiBi + d'ycCiCk 
= 0 , 

and the sum of correlation phases 7c is zero as the sum¬ 
mation terms cancel out each other. 
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